arXiv:1506.02383vl [astro-ph.EP] 8 Jun 2015 


Astronomy & Astrophysics manuscript no. hdl00546 pinilla © ESO 2015 

June 9, 2015 


Sequential planet formation in the HD 100546 protoplanetary disk? 

P. Pinilla\ T. Bimstiel^ and C. Walsh^ 

* Leiden Observatory, Leiden University, P.O. Box 9513, 2300RA Leiden, The Netherlands 
e-mail: pinilla@strw.leidenuniv.nl 

^ Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA 
Accepted 25 May 2015 


ABSTRACT 


Context. The disk around the Herbig Ae star, HD 100546, shows structures that suggest the presence of two companions in the disk 
at ^ 10 and ~ 70 AU. The outer companion seems to be in the act of formation. 

Aims. Our aims are to provide constraints on the age of the planets in HD 100546 and to explore the potential evidence for sequential 
planet formation in transition disks such as HD 100546. 

Methods. We compare the recent resolved continuum observations of the disk around HD 100546 with the results of dust evolution 
simulations using an analytical prescription for the shapes of gaps carved by massive planets. 

Results. An inner pressure bump must have been present since early in the disk lifetime to have good agreement between the dust 
evolution models and the continuum observations of HD 100546. This pressure bump may have resulted from the presence of a very 
massive planet (~ 20 Mj„p), which formed early in the inner disk (r ~10 AU). If only this single planet exists, the disk is likely to 
be old, comparable to the stellar age (~5-10 Myr). Another possible explanation is an additional massive planet in the outer disk 
(r ~ 70 AU): either a low-mass outer planet (< 5 Mj^p) injected at early times, or a higher mass outer planet (> 15 Mj^p) formed very 
recently, traps the right amount of dust in pressure bumps to reproduce the observations. In the latter case, the disk could be much 
younger (~ 3.0 Myr). 

Conclusions. In the case in which two massive companions are embedded in the disk around HD 100546, as suggested in the literature, 
the outer companion could be at least >2.5 Myr younger than the inner companion. 
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1. Introduction 

Transition disks display different interesting structures such as 
dust depleted cavities (e.g. Brown et al. 2009; Andrews et al. 
2011), azimuthal asymmetries (e.g Isella et al. 2013; Casassus 
et al. 2013; van der Marel et al. 2013; Perez et al. 2014), spi¬ 
ral arms (e.g. Fukagawa et al. 2004; Muto et al. 2012; Grady 
et al. 2013), and spatial segregation of small and large particles 
(e.g. Follette et al. 2013; Garufi et al. 2013; Zhang et al. 2014; 
Hashimoto et al. 2014). These structures suggest that these disks 
host a massive planet or multiple planets. Observations of planet 
candidates in transition disks (e.g Huelamo et al. 2011; Kraus & 
Ireland 2012; Biller et al. 2012; Quanz et al. 2013) have further 
supported this idea. Nonetheless, other mechanisms such as pho¬ 
toevaporation may also play an important role and explain some 
of the observed structures (e.g. Alexander et al. 2006; Owen & 
Clarke 2012; Rosotti et al. 2013). 

Observations of the disk around the Herbig Ae star 
HD 100546 have indicated the presence of two potential com¬ 
panions. By modelling the [O I] 6300 A emission line obtained 
with VFT/UVES in HD 100546, Acke & van den Ancker (2006) 
suggested the presence of a massive planet (> 20 Mjup) at 
~ 6.5 AU distance from the star. Using VLTI/AMBER at H- 
and K-band, Tatulli et al. (201 1) proposed a less massive planet 
(~ 1 - 8 Mjup) in the inner disk. However, using MIDI/VLT ob¬ 
servations and by constraining the curvature of the disk wall of 
the inner cavity, Mulders et al. (2013) found a lower limit for 
the mass of the inner companion. Taking temperature changes 
when a planet opens a gap into account, Mulders et al. (2013) 
conclude that the mass of a potential planet located at ~ lOAU 


is 20 - 30 Mjup. On the other hand, high-contrast imaging with 
VFT/NACO shows signatures of a massive planet (~ 15 Mjup) 
at ~ 70AU, which may be its formation stage (Quanz et al. 
2013, 2014). Currie et al. (2014) reported Gemini/NICI thermal 
infrared data and detected the outer protoplanet at the same lo¬ 
cation and brightness as that found by Quanz et al. (2013). This 
emission seems to have an extended structure which may come 
from a circumplanetary disk. In addition, spiral arms have been 
observed in scattered light images (e.g. Grady et al. 2001; Ardila 
et al. 2007; Boccaletti et al. 2013; Avenhaus et al. 2014), which 
may be related to the presence of these planets. 

Previous observations of protoplanetary disks have revealed 
that the disk radial extent can be much larger for the molecular 
gas than for the millimetre-sized particles. This is the case for 
TW Hya, whose CO emission extends up to ~215AU, while the 
mm grains extend to only ~60AU from the star (Andrews et al. 
2012). Bimstiel & Andrews (2014) suggest that radial drift may 
be responsible for the different radial extents of the gas and large 
dust grains. However, the disk around HD 100546 appears to be 
a special case: the radial extent of the gas is ~ 400AU, while 
most of the emission at millimetre wavelengths comes from a 
narrow ring concentrated at ~ 26AU, with a width of ~21 AU. 
In addition, a much fainter ring of emission (a factor of ~100 
lower) comes from the outer disk, which is centred at ~ 190 AU, 
with a width of ~ 75 AU (Walsh et al. 2014). This double-ring 
emission is consistent with the two-planet scenario (Walsh et al. 
2014). Ring-like emission in transition disks can be explained by 
the dust evolution, which occurs when a single massive planet or 
multiple planets interact with the disk (e.g. Pinilla et al. 2012a, 
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2015; Zhu & Stone 2014). An alternative explanation for multi¬ 
ple rings is magneto-rotational instability in the outer regions of 
disks (e.g. Pinilla et al. 2012b; Flock et al. 2014). In protoplane¬ 
tary disks, particles migrate inward because of the sub-Keplerian 
motion of the gas (e.g. Weidenschilling 1977; Nakagawa et al. 
1981; Tanga et al. 1996; Birnstiel et al. 2010). When a massive 
planet opens a gap, particles stop their inward migration due to 
the positive pressure gradient at the outer edge of that gap. Dust 
grains accumulate and grow in this preferential region known as 
a pressure trap, whose location and structure depend on disk vis¬ 
cosity, the mass and location of the planet (e.g. Rice et al. 2006). 

In this work, we investigate the influence of the two sug¬ 
gested planets in HD 100546 on the dust distribution in the disk, 
and compare the resulting predicted continuum emission with 
the most recent mm observations. By computing dust evolution 
models for different parameters, we aim to put constraints on 
different planet properties, and to address the following ques¬ 
tions: can dust evolution, without any planet(s) in the disk, ex¬ 
plain the two-component emission observed with ALMA? What 
happens when a single inner companion is assumed? What if two 
companions are assumed? To answer these questions, we cover 
a large parameter space, which includes the mass, location, and 
age of the planets, and disk viscosity. '. 

For HD 100546, proper hydrodynamical simulations for the 
planet-disk interaction are highly computationally demanding 
because of the large radial separation between the two planet 
candidates. Since we aim to study several cases with different 
planet properties, we instead use analytical solutions to model 
the shapes of gaps in the disk. We use the Crida et al. (2006) 
prescription for the width, and the Fung et al. (2014) prescrip¬ 
tion for the depth. In Sect. 2, we explain the assumptions for 
the carved gaps, the potential caveats, and the connection with 
dust evolution models. In Sect. 3 we present the results of the 
dust evolution models, the computed visibilities at different 
wavelengths (A - [0.87,1.0,3.0,7.0] mm), and the com¬ 
parison with ALMA and ATCA observations of HD 100546. 
Sections 4 and 5 are the discussion and conclusions of this work. 


2. Models and set-up 

In this section, the analytical approximation for the shape of the 
gaps carved by planets is described, together with the dust evo¬ 
lution models, and the computation of the visibilities. To test this 
analytical approach and the validity of the resulting dust density 
distributions, different benchmark cases are considered. 

2.1. Carved gaps 

Because of the large separation between the potential planets 
in HD 100546 (~ lOAU and ~ 70AU), hydrodynamical sim¬ 
ulations are computationally expensive. High grid resolution is 
needed close to the planet locations, and the inner and outer ra¬ 
dial boundaries should be far enough from the planet positions 
to avoid unphysical results. Because our interest is focussed on 
the influence of the gas surface density carved by planets on 
the radial dust evolution, we use an approximation of the radial 
shape of the gaps in protoplanetary disks rather than full hy¬ 
drodynamical simulations. We assume the analytical results pre¬ 
sented in Crida et al. (2006) and Fung et al. (2014). Crida et al. 
(2006) proposed an equilibrium profile considering the viscous 

* By age, we mean the time since the injection of the planet into the 
simulations 


torque (fy), gravitational torque (tg), and the torque removed by 
pressure supported waves or pressure torque (tp). This analyti¬ 
cal profile therefore satisfies that ty + tg + tp - 0. In an cr-disk 
(Shakura & Sunyaev 1973), where disk viscosity is parametrised 
as V = ac^ / Q., the equilibrium solution for the gas surface den¬ 
sity (Z) when a planet located at is interacting with the disk 
(see Crida et al. 2006, Eqs. 11, 13, and 14) is 


tg - |ac^ 


( 1 ) 


with rn being the Hill radius of the planet (rn - rp{qf3y^^, 
where q is the planet-stellar mass ratio Mp/M*, and the planet 
orbital radii), Qp the angular orbital velocity of the planet, ///r 
the aspect ratio, Cs the sound speed, and a" a dimensionless func¬ 
tion given by 



(r-rp) 
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+ 200 
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rn 

rn 



( 2 ) 


This dimensionless function a" is an ansatz from 2D verti¬ 
cally isothermal simulations. When Eq. (1) is solved, a boundary 
condition needs to be imposed, which in our case is assumed to 
be the unperturbed density. Because the equilibrium solution of 
Eq. (1) assumes the gravitational and pressure torques to be null 
close to the planet (-2r// < r - Vp < 2rpi), the depth of the gap 
is not perfectly constrained. 

An analytical solution for the depth of the gap close to the 
planet is difficult to calculate because of the strong tidal grav¬ 
itational field of the planet in this region. Using two indepen¬ 
dent codes (PEnGUIn and ZEUS90), Fung et al. (2014) provided 
an empirical relation of the depth of the gap carved by a non¬ 
migrating giant planet in a locally isothermal disk. They consider 
a large parameter space for the planet-stellar mass ratio q, vis¬ 
cosity a, and aspect ratio ///r, finding very similar results with 
both codes. 

We imitate the gas gaps carved by massive planets by solving 
Eq (1) and correcting the depth of the gap with the empirical 
scalings (equations 12 and 14 from Eung et al. 2014), between 
-Ivh < r - Kp < Ikh- Both the width and the depth of the gaps 
are essential for the final dust distributions in disks. The width 
determines the location of the pressure maximum at the outer 
edge of the gap and therefore the location of the peak of the 
millimetre emission (e.g. Pinilla et al. 2012a). The depth is also 
important because of the possible dust filtration (e.g Rice et al. 
2006; Zhu et al. 2012). 

To compare the results from these approximations, Eig. 1 
displays the comparison of the gas surface density between hy¬ 
drodynamical simulations done with FARGO (Masset 2000), and 
the gap shape obtained using the solution by Crida et al. (2006) 
and corrected by the empirical relations from Eung et al. (2014). 
For this comparison, it is considered that a massive planet, with 
q - I X 10“^ (left panel) and q - \ x 10“^ (right panel), opens 
a gap in a flared disk, i.e. H/r - liQrf, with /iq = 0.05 being the 
aspect ratio at the position of the planet, and / a flaring index 
equal to 0.25. The disk viscosity is assumed to be a = 10^^. 
From the hydrodynamical simulations, the gas surface density is 
azimuthally averaged after the disk reaches a steady-state, 1000 
and 2000 planet orbits for q' = lx 10“^ and q - \ x 10“^ re¬ 
spectively. As shown in Fig. 1, this approximation for the shape 
of the gaps is very good for the case of q' = 1 x 10“^. In the case 
of q - lx 10“^, the width is under-predicted compared to the 
azimuthally-averaged profile from the hydrodynamical simula¬ 
tions. When a very massive planet q > 5x 10“^ interacts with the 
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r{rp) r{rp) 

Fig. 1. Comparison between the analytical gap prescriptions and hydro-simulations when a massive planet opens a gap. The disk is 
assumed to be a flared disk, i.e. Hjr - h^rf , with /iq = 0.05 the aspect ratio at the position of the planet, and / flaring index equal 
to 0.25. The planet-to-star mass ratio k q - lx 10“^ (left panel) and q = 1 x 10“^ (right panel), and the planet position is - 1. 
The disk viscosity is assumed to be a = 10“^ for both cases. 


disk, the formed gaps become eccentric (Kley & Dirksen 2006), 
making the radial gap profile wider and less steep. The location 
of the pressure bump at the outer edge of the gap may change 
from 7-8 Rh to 9 - 10 Rh for a non-eccentric to an eccentric 
gap. Thus, in the case of transition disks, to reproduce the peak of 
the millimetre emission at a certain distance from the star using 
this analytical approach, we give an approximate value for the 
planet location rather than a specific position when q > 5x 10“^. 
Since the match between the hydrodynamical simulations and 
the analytical approach only works for the gap shape, and po¬ 
tential differences become significant far from the gap location 
(Fig. 1), we assume the gas surface density to follow the un¬ 
perturbed density far from the location of the planet. Additional 
azimuthal features that may exist when a massive planet inter¬ 
acts with the disk, such as vortices or spiral arms (e.g. Kley 
& Dirksen 2006; Ataiee et al. 2013; Fung et al. 2014; Zhu & 
Stone 2014; Juhasz et al. 2014), are not considered in this work. 
Instead, we focus on the radial distribution of particles. 


2.2. Dust evolution and radial gas velocity 


To model the dust evolution, we use the formulation of Birnstiel 
et al. (2010). This model solves the advection-diffusion differen¬ 
tial equation for the dust surface density and simultaneously 
simulates the growth, fragmentation, and erosion of dust grains 
by considering collisions of particles. For the relative velocities 
between particles. Brownian motion, turbulent velocities, set¬ 
tling to the midplane, and radial dust velocities are taken into 
account. The detailed explanation of this dust evolution model is 
in Birnstiel et al. (2010). This model has been used extensively to 
investigate dust distributions in different types of disks, includ¬ 
ing comparisons with observations (e.g. Birnstiel et al. 2013; 
van der Marel et al. 2013; Pinilla et al. 2013; de Juan Ovelar 
et al. 2013). We investigate the evolution of 180 species of dust 
grains, defined by size from 1 pm to 200 cm. The radial velocity 
of dust depends on the coupling of the dust particles to the gas. 
The Stokes number (St) quantifies this coupling, which in the 
Epstein regime, where the mean free path of the gas molecules 


dmfp is higher than the size of particles a (Tmfp > 4/9a), is de¬ 
fined as 


St = 


aps TT 

S 2’ 


(3) 


where ps is the volume density of a dust grain of size a. In terms 
of the Stokes number, the total radial dust velocity is 


Wg 1 drP 

1+St2 sr'-HStPgfJ’ 


(4) 


with pg being the total gas surface density. The first term of 
Eq. (4) depends on the radial gas velocity which in a disk 
evolving by viscous accretion is 


•"r.g 


3 d 


(Sv ^^r). 


(5) 


When the dust evolution models are combined with hy¬ 
drodynamical simulations, the gas surface density and gas ra¬ 
dial velocities can be directly taken as initial conditions for the 
dust evolution models, once the disk has reached a steady-state. 
Alternatively, the gas surface density can be based on the hydro- 
dynamical simulations and the gas radial velocity from Eq. (5) 
as in Pinilla et al. (2012a). When a non-migrating planet is in 
the disk, Y is calculated with the analytical solution described in 
Sect. 2.1 and using this E -profile, the gas radial velocities are 
calculated assuming Eq. (5). To see the potential differences of 
assuming the gas radial velocities from hydrodynamical simula¬ 
tions or viscous accretion, Eig. 2 compares the gas radial veloc¬ 
ities obtained by hydrodynamical simulations using FARGO for 
q - 1 X 10~^, Hjr - hor^ with ho = 0.05 and / = 0.25, and 
a - 10“^ (azimuthally averaged and averaging over the last 100 
orbits of evolution i.e orbits 1000-1100), when Ui d is calculated 
from viscous accretion (Eq. 5) and assuming Y from FARGO; and 
when the gas surface density is calculated analytically (Sect. 2.1) 
and Vr,d fromEq. (5). In addition the dust density distributions af¬ 
ter ~ 1 Myr of evolution are displayed for each case. The main 
differences of the gas radial velocities are close to the planet 
where strong perturbation waves due to the planet create a wig¬ 
gle profile. However, the dust density distribution is similar for 
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Fig. 2. Left panel: comparison of the gas radial velocity obtained by hydrodynamical simulations using FARGO for q - \ x 10“^, 
Hjr = h()r^ with ho - 0.05 and / = 0.25, and a = 10“^, when ^ is calculated from viscous accretion (Eq. 5), assuming S by 
azimuthally averaging the gas surface density from FARGO simulations, and when the gas surface density is calculated analytically 
(Sect. 2.1) and from Eq. (5). Other panels: dust density distribution at the same time of evolution (~ IMyr) for the same 
assumptions. White line corresponds to St = 1 (Eq. 3), which is proportional to the gas surface density. 


each case, because radial drift and diffusion are more important 
than the drag term in Eq. (4). Hence, we have also proved the 
reliability of this analytical approach compared to the hydrody¬ 
namical simulations to obtain proper radial density distribution 
of dust particles. 

2.3. Visibilities at mm-wavelenghts 

To compare the results from the dust evolution models with mil¬ 
limetre observations, we calculate the real part of the visibilities 
in the Mv-plane: VReai(r‘«v) is given by (Berger & Segransan 2007) 

/^OC 

VReai(^uv) = 27t j I(r)Jo(27rruv)rdr, (6) 

Jo 

where Jq is the zeroth-order Bessel function of the first kind and 
I{r) is the radial-dependent emergent intensity, which is directly 
calculated using the vertically integrated dust density distribu¬ 
tion cr(r, a) from the dust evolution results at a specific time of 
evolution. Thus, cr(r, a) is different at each location for all the 
dust particle sizes (180 species) assumed in this work. For a 
given wavelength (A), the intensity is given by, 

/(r) = B,(r(r))[l-exp-^^W], (7) 

with B;[(T(r)) being the Planck function at the temperature T{r) 
and Tx the optical depth, which is computed as 

= 2 ^, , 8 ) 

COS I 

where the opacities at a particular wavelength k,\ are calcu¬ 
lated for each grain size, assuming Mie theory and a mix of 
magnesium-iron silicates (e.g. Dorschner et al. 1995). The op¬ 
tical constants are taken from the Jena database 

2.4. Set-up 

We assume the disk mass to be 0.05 Mq, consistent with the 
dust mass from Mulders et al. (2013), and a canonical dust-to¬ 
gas mass ratio of 100; however, disk masses are very uncer¬ 
tain. From optically thin mm emission, the disk mass can be 

^ http://www.astro.uni- jena.de/Laboratory/Database/databases.html, 
with a specific silicate composition of: 10% MgFeSi 04 , 28% MgSiOs, 
31% Mg 2 Si 04 , 1% NaAlSizOfi. 


Table 1. Stellar, disk and planets parameters 


Parameter 

Symbol /units 

Value 

Stellar mass 

44* [Mel 

2.4 

Disk mass 

44disk[44o] 

0.05 

Inner disk radius 

fJAU] 

2.0 

Outer disk radius 

ro.dAU] 

400 

Viscosity 

a 

[2x 10-^5x 10-3] 

Inner planet mass 

44pi[44jiip] 

[1,5,10,20,30] 

Outer planet mass 

^/j2[^Jup] 

[5,10,15,20] 

Inner planet position 

rp,[AU] 

~ 10 

Outer planet position 

VlAU] 

~70 

Distance to the disk 

d[pc] 

103 

Disk inclination 

in 

45 

Fragmentation velocity 

Vf[m s“'] 

10 

Volume density of dust 

Pdg cm^^] 

1.2 


estimated assuming dust opacities and a dust-to-gas mass ratio 
(e.g. Andrews & Williams 2005). Nonetheless, this calculation 
is unreliable since the dust opacity depends on composition and 
shape of the dust (e.g. Min 2005; Demyk et al. 2013), and the gas 
and dust are not necessarily co-spatial (e.g. Birnstiel & Andrews 
2014). The gas mass can also be calculated from observations 
of CO and its isotopologues (e.g. Williams & Best 2014); how¬ 
ever, this estimation can depend on the chemical disk evolution, 
such as isotope selective processes (Miotello et al. 2014). The 
disk mass assumed for HD 100546 is broadly consistent with 
the mass estimate from e.g. Henning et al. (1998), but it remains 
a very uncertain parameter and thus we keep it fixed in this work. 

We fix the gas disk extent to [2 - 400] AU, in agreement with 
the CO J = 3-2 emission from Walsh et al. (2014). The initial 
gas surface density is a power law, such that S(r) = 2o(^/^/)i)~'- 
We use a temperature profile specific for HD 100546 and de¬ 
rived by Bruderer et al. (2012), who constrained the gas temper¬ 
ature via detailed modelling of the observed low-J CO line emis¬ 
sion from single dish observations with the Atacama Pathfinder 
Experiment (APEX) (Panic et al. 2010), mi d/high-J CO lines 
observed with PACS-Herschel (Sturm et al. 2010), and contin¬ 
uum emission. For the dust evolution models, we impose the 
dust temperature as the gas temperature in the midplane, where 
Tgas - Tdust (with values of T 20 AV - lOOK, TdooAu - 25K). This 
temperature is an upper limit, since it is obtained from obser¬ 
vations of the warm gas in the disk atmosphere. However, CO 
channel maps from ALMA observations (Walsh et al. 2014) do 
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Fig. 3. Dust density distribution at different times of evolution 
{top panel) and real part of the visibilities at 870/im {bottom 
panel) at different times of dust evolution, when no planet is em¬ 
bedded in the disk. ALMA Cycle 0 data at this wavelength are 
over-plotted (data from Walsh et al. 2014). In the upper panel, 
the white and blue line correspond to Adrift (Eq. 10), and a^ag 
(Eq.9), respectively. 


not show the double lobe signature, which indicates the presence 
of a CO freeze-out zone (see e.g. de Gregorio-Monsalvo et al. 
2013). Hence, HD 100546 is a warm disk (Lniidpiane ^ 20 - 25K 
everywhere). Thus, we do not expect that the disk temperature 
considerably differs from the adopted temperature. This tem¬ 
perature is used for the dust evolution models and assumed to 
be the temperature of the vertically integrated dust density dis¬ 
tribution for the calculation of the intensity and visibilities at 
mm-wavelengths, which are mostly sensitive to the distribution 
of large grains. Because large (mm-sized) grains efhciently set¬ 
tle towards the midplane, the surface density of large grains in 
the disk atmosphere (where the temperature is higher) is neg¬ 
ligible (Dullemond & Dominik 2005). Therefore, adopting the 
midplane temperature only in the calculation of the emergent 
continuum emission is an appropriate assumption. 

Based on that temperature, the aspect ratio is calculated as 
CsIvk - Hjr - hor-f, obtaining Iiq - 0.045 and / = 0.33 at the 
position where the inner planet is located (10 AU). This aspect 
ratio is directly used for obtaining the analytical shape of carved 
gaps. The disk viscosity is assumed to be v = aCj/fl, with a - 


[2x10^^, 5x10“^] as in Mulders et al. (2013). For the visibilities, 
the disk is taken to be at a distance of 103 pc with an inclination 
(/) of 45° (van den Ancker et al. 1997), and a position angle 
(P.A.) of 146° (east from north). 

We assume a large range of planet masses according to 
previous studies (Acke & van den Ancker 2006; Tatulli et al. 
2011; Mulders et al. 2013; Quanz et al. 2013; Currie et al. 
2014). For an inner companion located ~ 10 AU, we as¬ 
sume the planet-stellar mass ratio to be ~ [4.2 x 10“^, 2.1 x 
10“^,4.2 X 10“^, 8.3 X 10“^, 1.3 x 10“^], which corresponds to 
Mpi - [1,5,10,20,30] Mjup planets around a 2.4 Mq star. For 
the outer planet located at r,,? ~ 70 AU, we consider M „2 = 
[5,10,15,20]Mj,p. 

Finally, for the dust evolution models, we assume that the 
velocity threshold above which particles fragment, the so-called 
fragmentation velocity, is Vf = 10 ms“'. For the initial dust 
density distribution, we assume all particles to be 1 yum-sized. 
The dust particles are considered to have a volume density of 
Ps - 1.2g cm“^, according to the averaged values of the volume 
density for silicates (e.g. Blum & Wurm 2008). All parameters 
are summarised in Table 1 . 


3. Results 

In this section, we present the results for the case of no planet, 
and for the cases where either one or two planets are embedded 
in the disk. For the latter case, we consider the two scenarios 
where the second (outer) planet is injected either at the same 
time as the inner planet, or at a later time. 

For the models where the planets are injected at later times, 
we do not assume an actual model to introduce the planet and 
let it grow. Our main aim is to study the final dust density dis¬ 
tributions assuming that the planets are already formed and em¬ 
bedded in the disk. Massive planets such as those assumed in 
this work are expected to have slow type II migration, in which 
case migration timescales follow the viscous diffusion time of 
the gas, as do the pressure bumps and the dust. Thus, potential 
dust traps are likely to follow the migration of the planets, keep¬ 
ing qualitatively similar results. 


3.1. No planets in the disk 

Figure 3 illustrates the dust density distribution at different times 
of evolution, and the corresponding real part of the visibilities 
at 870jum, when no planet is embedded in the disk. Moreover, 
ALMA Cycle 0 data at this wavelength are also plotted for com¬ 
parison. Note that the errors are also plotted and they are smaller 
than the point size (because of the very high signal-to-noise ra¬ 
tio of the ALMA data, Walsh et al. 2014). For this simulation 
a - 2x10“^ is assumed. At early times (f < 1 Myr), the dust par¬ 
ticles grow to the maximum grain size before particles fragment. 
When fragmentation is mainly because of turbulent relative ve¬ 
locities, Afrag is (Birnstiel et al. 2012), 


- 


(9) 


At later times, the dust surface density decreases, so the drift 
barrier moves to smaller sizes and possibly below the fragmen¬ 
tation size (Eq. 9). In such cases, Aduft is (Birnstiel et al. 2012) 
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Fig. 4. Dust density distribution after 1.5 Myr of evolution and real part of the visibilities at 870/tm at different times of dust 
evolution, when a single planet interacts with the disk for two different values of disk viscosity, a = 2 x 10~^ (left column) and 
a = 5x 10^^ (right column). From top to bottom, the masses of the planet are; 1, 5, and lOMjup. ALMA Cycle 0 data are over-plotted 
in Fig. 3. In the dust density distribution plots, the white and blue lines correspond to St = 1 (Eq. 3), and afmg (Eq. 9) respectively. 


Before particles start to drift, millimetre particles are dis¬ 
tributed in the entire disk (Eig. 3), creating a visibility profile in 
disagreement with observations, which suggest a more concen¬ 
trated ring centred at ~ 26AU (~ 290 kd) (Walsh et al. 2014). 
Once the particles grow to sizes for which radial drift dominates 
(f > 1 Myr), the disk is quickly depleted of millimetre dust par¬ 
ticles. This is reflected in the profile of the real part of the visi¬ 
bilities where the flux at 870jum drastically reduces after 1 Myr 
of evolution, becoming already very low at 2 Myr. 

As a consequence, any mechanism that helps to reduce the 
rapid inward drift is needed to explain the millimetre emis¬ 
sion of this disk. We assume that a pressure trap is formed at 
the outer edge of a gap carved by a massive planet, motivated 
also by the companion candidates suggest in the literature (e.g. 
Acke & van den Ancker 2006; Tatulli et al. 2011; Brittain et al. 
2014). Particles concentrate in the regions of highest pressure 
(e.g. Weidenschilling 1977; Brauer et al. 2008), and therefore 
the assumed pressure bump can help to reduce the radial drift. 

3.2. A planet is embedded in the inner disk 

Assuming the analytical profiles for gaps carved by massive 
planets described in Sect. 2.1, we consider an inner planet of 
different masses Mpi — [1,5,10,20,30] Mjup at ~ lOAU to 
investigate the trapping efficiency and the comparison with ob¬ 
servations. Since we neglect the possible effect of an eccentric 
gap in our simulations, we give an approximated value for the 
location of the planets (Sect 2.1). 

Eigure 4 shows the dust density distribution after 1 Myr of 
evolution, and the corresponding real part of the visibilities at 
870//m for a = [2 x 10“^, 5 x 10“^], when a 1, 5 or 10 Mjup 


planet is embedded in the inner part of the disk. The results of 
these simulations suggest that a planet whose mass is lower than 
10 Mjup is not massive enough to result in trapping of millimetre¬ 
sized dust particles, sufficient to create a visibility curve that is 
in agreement with ALMA Cycle 0 observations. In any of these 
cases, the flux is already under-predicted after 1.5 Myr of evolu¬ 
tion. Pinilla et al. (2012a) show that when a planet of ^ = 10“^, 
i.e. a 1 Mjup around a solar-type star, carves a gap in the disk, 
particles can be trapped at the outer edge when a = 10“^. The 
reason why trapping is more difficult in this case is because of 
the larger difference between the maximum grain size or afmg 
(Eq. 9), and the particle size in which dust feels the highest ra¬ 
dial drift i.e. St = 1. Particles with St ~ 1 are the easiest to trap, 
however a^ag is in this case around one order of magnitude lower 
than the grain size that corresponds to St = 1. These smaller par¬ 
ticles are more difficult to trap because the drift is inefficient for 
particles with St < a. Eor these models, the maximum grain size 
is lower because of the higher temperatures and a-viscosity con¬ 
sidered here compared with those in Pinilla et al. (2012a). The 
visibility shape does not change significantly for the two val¬ 
ues of the viscosity. For this reason, for the following results, we 
only focus on the cases with a - 2x10^^. Increasing a decreases 
the maximum grain size even further (Eq. 9), making the particle 
trapping more difficult. Assuming a lower viscosity or a higher 
fragmentation velocity can help to trap the particles (Pinilla et 
al. 2015). 

When the mass of the inner planet is increased, the trapping 
becomes more efficient. Eigure 5 shows the dust density distribu¬ 
tion at different times of evolution and the real part of visibilities 
at T = [0.87,1.0,3.0,7.0]mm, when a 20 and 30 Mjup planet is 
embedded in the inner part of the disk. ALMA Cycle 0 data from 
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Fig. 5. Dust density distribution at different times of evolution {upper row) and the corresponding real part of the visibilities at 
A — [0.87,1.0,3.0,7.0] mm {lower row) at different times of dust evolution, when a 20 Mjup planet {left column) and 30 Mjup planet 
{right column) is embedded in the inner disk. Data from ALMA Cycle 0 from Walsh et al. (2014) and ATCA from Wright et al. 
(2015) are over-plotted. In the upper panels, the white and blue line correspond to St = 1 (Eq. 3), and (Eq.9) respectively. 


Walsh et al. (2014) and ATCA data from Wright et al. (2015) are 
over-plotted for comparison. Eor the case of a 20 Mjup planet, 
the visibilities better ht the data; nonetheless, to obtain the best 
ht with the data, the dust needs to evolve for long timescales 
(> 5 Myr) to grow, drift, and concentrate the millimetre grains in 
a narrow ring. Already at 5 Myr, most of the millimetre particles 
are concentrated in a narrow region; however, the large amount 
of these grains leads to a over-prediction of the fluxes. The con¬ 
tinuous fragmentation that occurs in the pressure bump, together 
with some small-sized dust crossing the gap, reduce the amount 
of millimetre grains, enabling a better ht to the data at 10 Myr. 
The total hux is under-predicted at 7mm because the maximum 
grain size in the trap is around 1mm. The null of the visibilities 
at ~ 290 k/1 is in good agreement with the observations, imply¬ 
ing that a 20 Mjup planet located in the inner part of the disk, can 
create a pressure maximum at ~26AU, where the centre of the 
inner narrow ring was observed with ALMA Cycle 0. 

If the planet is located at the same location, and its mass is 
increased to 30 Mjup, the resulting carved gap is wider, moving 
the location of the pressure maximum outwards and therefore the 
peak of the millimetre emission. Eor this reason, the null of the 
visibilities is at shorter baselines (~ 245 kT, Eig. 5). In addition. 


the total flux is higher because the gap is also deeper, hltering 
more dust, and as a result, a larger number of millimetre parti¬ 
cles remain in the trap. However, the fluxes at 7mm are lower 
than in the case of a 20 Mjup planet. This is because the maxi¬ 
mum grain size has decreased slightly in the trap as the pressure 
maximum is further away, where the gas surface density is also 
lower, decreasing (Eq. 9) within the bumps. An important 
remark is that at 7 mm the emission likely has a > 10% con¬ 
tamination from free-free emission (Wright et al. 2015), which 
is neglected in our calculations. Hence, we expect the models to 
underestimate the 7 mm flux. 


Erom the simulations with one single planet in the inner 
disk, the best model that reproduces the observations at different 
wavelength is for a 20 Mj^p planet and a - 2 x 10“^, consistent 
with the results from Mulders et al. (2013), who found similar 
results from modelling mid-infrared data. Eor these parameters, 
the best fit is for an old disk (~ 5 - 10 Myr). Nevertheless, with 
one single planet, the outer ring of mm emission observed with 
ALMA remains unexplained (Walsh et al. 2014) . 
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Fig. 6. Upper panels: dust density distribution at 1.5 Myr of evolution when two planets are injected in the disk simultaneously. The 
mass of the inner planet is 20 Mjup and it is located at ~ 10 - 15AU, while the outer planet at rp 2 >70 AU is assumed to be 
^pi - [5,10,15,20] Mjup from left to right panels respectively. The white line is St = 1, and the blue line is ctfrag. The corresponding 
real part of the visibilities for each Mp 2 at T = 0.87mm {middle panels) and A - 7.0mm {lower panels) are displayed for different 
times of dust evolution, for each Mp 2 . ALMA Cycle 0 and ATCA data are over-plotted for comparison. 


3.3. Two planets are embedded in the disk 

In this section, the results with two planets embedded in the disk 
are presented. We use the analytical shapes described in Sect. 2.1 
for both gaps. For all cases, the gas surface density is always 
scaled such that the disk mass is 0.05 Mq. For the planet em¬ 
bedded in the inner disk, we consider the best single-planet fit of 
Mpi - 20 Mjup, and a disk viscosity of a = 2x 10“^. For the outer 
planet, it is assumed Mp 2 - [5,10,15,20] Mjup and rp 2 ~ 70AU. 
There are two sets of simulations: one where the two planets are 
assumed to interact with the disk from the same stage of evolu¬ 
tion, and a second, where the outer planet is assumed to interact 
with the disk from a later stage. 

3.3.1. Simultaneous injection of both planets 

Figure 6 illustrates the dust density distribution at 1.5 Myr of 
evolution when two planets are interacting in the disk from the 
same stage of evolution (f ~ 1000 yr). The inner planet is 
a 20 Mjup planet at ~ 10 AU, while the outer planet at 

rp 2 ~70 AU is assumed to be Mp 2 - [5,10,15,20] Mjup. In 
addition, the real part of the visibilities at T = [0.87,7.0]mm 
are displayed for each Mp 2 , and at different times of dust evo¬ 
lution. When a second planet is assumed, the gas surface den¬ 
sity sharply decreases with radius in the region between the two 
planets. As a consequence of the high negative pressure gradi¬ 
ent, the dust particles move quickly inward (timescales shorter 
than ~ 1 Myr) and stop their migration in the pressure maximum 


at the outer edge on the inner gap. This effect creates a narrow 
ring-like concentration of large particles centred at ~ 26AU in 
shorter timescales than in the case of a single inner planet. 

In this case, a second trap also exists at the outer edge of 
the gap carved by the planet at rp 2 ~70 AU. The amount of 
dust concentrated in this second trap depends on the mass of 
the outer planet. When the planet is more massive, dust filtration 
becomes more effective and the millimetre flux increases. This 
has a significant influence on the resulting real part of the visi¬ 
bilities. Because of the dust trapping in both pressure maxima, 
the visibility profiles show strong oscillations. For the case of 
Mp 2 = 5 Mjup, these oscillations are smoothed out late in the 
evolution (~10 Myr) because of the continued crossing of dust 
particles through the outer gap, which is filling the inner trap 
(middle and bottom left panels of Fig. 6). In this case, there is 
a similar contrast of the inner to the outer ring emission as that 
observed with ALMA Cycle 0 (of the order of ~100), generat¬ 
ing similar visibility profiles to those observed. Nevertheless, for 
the case of a higher mass planet {Mp 2 > 10 Mjup), the amount 
of dust concentrated in both traps is similar, and the undulating 
shape is more extreme and persists at late times, in contradic¬ 
tion to observations. The total flux at T = 0.87 mm is slightly 
lower than in the case of the 20 Mjup single planet; however, 
these fluxes can be higher if the disk mass or dust-to-gas ratio 
are assumed to be higher, or they can change by assuming dif¬ 
ferent dust composition (Min 2005). We keep these two values 
the same for all the simulations to have a clear distinction of the 
effect of the planet parameters on the final dust distributions. 
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Fig. 7. Upper panels: dust density distribution at 2.5 Myr of evolution when a 20 Mj^p planet is embedded in the inner disk (10- 
15AU) and a second planet (15 Mjup left panels, 20Mjup right panels) is injected in the outer region {rp 2 >70 AU) after 2 Myr of dust 
evolution. The white line is St = 1, and the blue line is flfrag. Bottom panels: real part of the visibilities at d = [0.87,1.0,3.0,7.0]mm 
at different times of dust evolution; ALMA Cycle 0 and ATCA data are over-plotted for comparison. 


From these two-planets simulations, in which both planets 
are considered to interact with the disk from early in the disk 
lifetime, we conclude that the outer planet should be a low-mass 
planet (< 5 Mjup) compared to the inner planet (20 Mjup), and 
the disk may be as old as in the single planet scenario (~5- 
10 Myr). 

3.3.2. Later injection of the outer planet 

High-contrast imaging of this disk with VLT/NACO shows sig¬ 
natures of a massive planet (~ 15 - 20 Mjup) in the outer region 
at ~ 70 AU (Quanz et al. 2013; Currie et al. 2014). However, 
in the previous section we showed that if the planet in the outer 
disk is as old as the one in the inner disk, the resulting visibility 
profiles are in disagreement with observations. To have a very 
massive planet in the outer disk and similar millimetre emission 
contrast between the inner and the outer ring (~ 100), the outer 
planet must be younger than the inner planet. In this scenario, at 
the time that the outer planet is injected in the disk, most of the 
dust has already moved towards the inner trap, decreasing the 
mass of dust available for trapping in the second pressure bump. 

Figure 7 shows the dust density distribution at 2.5 Myr of 
evolution when a 20 Mjup planet is embedded in the inner disk, 
and a second planet {Mp 2 - [15,20] Mjup) is injected in the outer 
region {rp 2 ~70 AU) after 2 Myr of dust evolution. In compar¬ 
ison with the corresponding cases of Fig. 6, it is important to 


notice that the outer ring is fainter in this case than when the 
two planets are assumed to interact with the disk over the same 
timescale. The oscillations on the visibility profiles are much 
smoother in these cases, and better fit the total flux. The null 
in the visibility profiles remains similar as in the single planet 
case (~ 290 kd). A 15 Mjup outer planet fits the visibility shapes 
better and total fluxes than a 20 Mjup planet, in agreement with 
the Gemini/NICI observations reported by Currie et al. (2014), 
who constrain the mass of the outer planet to be at most 15Mjup. 
The real part of the visibilities still have smooth oscillations; 
however, this can be further smoothed out by assuming than the 
planet is injected even later in the simulations. Figure 8 displays 
the case of a 20 Mjup planet embedded in the inner disk from 
early stages, and an outer 15 Mjup planet injected after 3 Myr 
of dust evolution. To reproduce the observations, the lower the 
mass of the planet, the earlier it needs to be introduced into the 
outer disk. Therefore, assuming the mass predicted by Quanz et 
al. (2013, 2014) and Currie et al. (2014) for the outer planet (and 
its possible circumplanetary disk), implies that this planet must 
be at least >2-3 Myr younger than the planet in the inner disk. 

4. Discussion 

We have shown that radial drift alone cannot reproduce the dust 
distribution observed for HD 100546 (Sect. 3.1). To obtain good 
agreement between the dust evolution models, which include 
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radial drift, with the millimetre observations, an inner pressure 
bump must exist since early in the disk lifetime (<0.1Myr). 

Observations of HD 100546 indicate the presence of at least 
two companions in this disk (see e.g. Acke & van den Ancker 
2006; Quanz et al. 2013). The mass of the inner companion 
has largely been discussed by several authors (e.g. Tatulli et 
al. 2011). If the required inner pressure bump is induced by a 
planet, we demonstrated that the required planet must be massive 
(~ 20 Mjup) for a disk viscosity of a = 2 x 10“^. This in agree¬ 
ment with Mulders et al. (2013) results, but is derived from dif¬ 
ferent observations (Sect 3.2). The trapping in HD 100546 disk 
is difficult because of the high disk temperature (as derived in 
Bruderer et al. 2012) and disk viscosity, which both increase the 
relative motion of the particles due to turbulence, and decrease 
the maximum grain size that particles can reach under these as¬ 
sumptions. This maximum grain size is one order of magnitude 
lower than the size corresponding to St = 1, which are the parti¬ 
cles that feel the highest radial drift and are the easiest to trap in 
a pressure bump. Assuming a higher viscosity would make the 
trapping even more difficult and a more massive inner planet, lo¬ 
cated closer to the star, would be needed to explain the millime¬ 
tre emission. In the case of a single planet embedded in the inner 
part of the disk, the disk is more likely to be old (~5-10 Myr). 

In addition to the inner planet, we also assumed an outer 
planet as suggested by various authors e.g. Boccaletti et al. 
(2013) and Quanz et al. (2013). When the two planets are as¬ 
sumed to interact with the disk from the same stage of evolu¬ 
tion, we demonstrated that the mass of the outer planet needs 
to be lower than the inner planet (< 5 Mjup), otherwise the re¬ 
sulting dust density distributions are in disagreement with mil¬ 
limetre observations. This is because, if a second massive planet 
is coeval with the inner planet, it traps inward drifting grains 
too effectively in the outer bump, and the contrast ratio be¬ 
tween the bumps becomes too strong compared with observa¬ 
tions. However, the mass suggested by various authors for the 
outer planet is higher than 5 Mjup. As in the single planet case, 
the models are more consistent with observations when the disk 
is old (~ 5-lOMyr). To test how sensitive the resulting timescales 
are to the assumed gas surface density profile, we considered the 
best-fit single planet model and adopted an exponentially tapered 
disk rather than a power law. For this case, mm-grains are con¬ 
centrated in a much narrower ring, shifting the null of the vis¬ 
ibilities to longer baselines, in disagreement with observations. 
In addition, at very early times, there is an oscillating behaviour 
(similar to the previous two-planet simulations). Hence, with this 
gas surface density profile, stronger oscillating behaviour is ex¬ 
pected for the visibilities for the case of two planets. Thus, a 
power law for the gas surface density gives us a better ht to the 
observations. 

If the mass of the outer planet is assumed to be large, as 
suggested by Quanz et al. (2013) (~ 15 Mjup), we demonstrated 
that the outer planet must be at least ~ 2 - 3 Myr younger than 
the inner planet, favouring a signihcantly younger outer planet 
observed in the act of formation. In this case, the disk can be 
much younger than when the disk only hosts a single planet. 
The required time to be in agreement with observations is a few 
tenths of a Myr after the outer planet is injected. Figure 9 sum¬ 
marises our main hndings and compares the cases of no planet, a 
single planet, two planets injected simultaneously, and two plan¬ 
ets where the outer planet is injected after 3 Myr of evolution. 
One of the uncertainties of these predictions is the disk mid¬ 
plane temperature, which may be slightly lower than obtained by 
Bruderer et al. (2012) (Sect. 2.4). The temperature prohle used is 
an upper limit; however, ALMA Cycle 0 CO J=3-2 observations 



Deprojected baseline (kA) Deprojected baseline (kA) 

Fig. 8. Real part of the visibilities at T = [0.87,1.0,3.0,7.0]mm 
at different times of dust evolution when a 20 Mjup planet is em¬ 
bedded in the inner disk (10-15AU) from early stages, and a sec¬ 
ond planet (15 Mjup) is injected after 3 Myr of dust evolution. 
ALMA Cycle 0 and ATCA data are over-plotted for compari¬ 
son. 
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Fig. 9. Real part of the visibilities at 870;um. Case 1: no planet 
embedded in the disk and f = 2 Myr. Case 2: a single planet em¬ 
bedded in the disk of 20Mji,p mass at ~ 10 AU and f = 10 Myr. 
Case 3: two planets simultaneously interacting with the disk, 
20Mjup at ~ 10 AU and 5Mjup at ~ 70 AU, and f = 10 Myr. 
Case 4: two planets in the disk, but the outer planet is injected 
t - 3Myr after the inner planet, 20Mjup at ~ 10 AU and 15Mjup 
at ~ 70 AU, and t - 3.5 Myr. ALMA Cycle 0 is over-plotted for 
comparison. 


show no signs of CO freeze-out, i.e. we also have a lower limit 
of ~20 K. A slightly lower midplane temperature decreases the 
dust drift velocities (Eq 4), and increase the maximum grain size 
(Eqs. 10 and 9). Observations with higher angular resolution and 
sensitivity of optically thin emission from multiple transitions of 
CO isotopologues are needed to give better constraints on the 
disk midplane temperature and confirm the current predictions. 
However, we expect that the overall trends remain similar with a 
slight lowering of the dust temperature. 

Moreover, Eig. 10 shows the synthetic ALMA Cycle 0 con¬ 
tour lines image at 870jum (imaged using identical (m, v) coor¬ 
dinates as the ALMA observations), for the cases of one single 
planet and an old disk, and two massive planets and a younger 
disk. The latter scenario reproduces a more extended emission 
than with a single planet, but it is not as extended as observed. 
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Fig. 10. ALMA Cycle 0 continuum contours at 870/tm (red 
dashed lines) and the contours corresponding to the models (blue 
solid lines) when an inner 20 Mjup planet alone (top panel) is 
embedded in the disk, and with two planets in the disk (bottom 
panel), where the outer planet is injected after 3 Myr of evolu¬ 
tion. Contours are every 3, 10, 30, 100, 300, and 1000 times the 
rms (0.5 mJy beam"'). 


This discrepancy is because the width of the carved gaps (thus 
the location of the pressure maximum at the outer edge of the 
gap) is underestimated in our models compared to proper hydro- 
dynamical simulations (Sect 2.1). The location of the pressure 
maximum of a gap carved by a 15 Mjup planet at 70 AU is ex¬ 
pected to be ~ 190 AU i.e. at 9 - lOr/^ from the planet location, 
where the peak of the outer ring emission is observed. Figure 11 
illustrates the synthetic images at 870 pm convolved with a 
0.1" X 0.1" beam from the same models used for Fig 10, show¬ 
ing that with high angular resolution observations, two rings can 
be resolved in the two-planet scenario. 
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Fig. 11. Synthetic images at 870 pm convolved with a O.F'xO.l" 
beam from the models of a single planet (20 Mjup at ~10 AU) 
and two planets (20 Mjup at ~10 AU and 15 Mjup at ~70 AU), 
injected 3 Myr after the inner planet. 
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The age of HD 100546 itself is very uncertain (~3-10Myr 
Acke & van den Ancker 2006) and a more precise measurement 
of the age of this star can give us hints as to the nature of the 
potential planets embedded in this disk. An observation to ul¬ 
timately test dust trapping in two pressure bumps, would be to 
measure spectral index variations inside and outside these pres¬ 
sure maxima locations (Pinilla et al. 2014a), which are testable 
with future ALMA capabilities. Figure 12 shows the expected 
radial variations of the spectral index, calculated between 1 and 
3 mm, when a single inner planet of 20 Mjup mass is embedded in 
the inner disk at ~10 AU, and when two planets are in the disk, 
but the outer planet (15 Mj^p) is injected 3Myr after the inner 
planet (20Mjup). At the location of the pressure trap, the spectral 
index is lower because of the accumulation of mm-grains and 
vice versa. 

We also used ATCA observations to compare the visibility 
profiles at 3 and 7 mm with the models. The models with two 
planets predict slightly better the total flux, and estimate well 
the observed null of the real part of the visibilities at these wave¬ 
lengths. The null at 7 mm is at a moderately longer kA than at 
shorter wavelengths. This happens because larger grains are ex¬ 
pected to be more concentrated. Nonetheless, for all cases, the 
total flux at 7mm is under-predicted, which partially comes from 
the neglect of free-free emission in the models, and because of 
the small maximum grain size (~ 1mm) under the assumed tem¬ 
perature and disk viscosity in the outer disk. 


5. Conclusion 

To obtain good agreement between dust evolution models and 
millimetre observations of the disk around HD 100546, an inner 
pressure trap must exist from early in the disk lifetime (<1 Myr). 
If this pressure bump is formed because of a planet embedded in 
the inner disk (r ~10 AU), the mass of the planet should be high 
(~ 20 Mjup) for a disk viscosity of a = 2 x 10"^, in agreement 
with Mulders et al. (2013), who found similar results from mod¬ 
elling mid-infrared data. In the case in which this is the only 
planet embedded in the disk, the disk is more likely to be old 
(~5-10 Myr). In this case the outer ring of emission observed 
with ALMA remains unexplained (Fig 10). If an outer planet is 
also embedded in the disk, and it is as massive as suggested by 
Quanz et al. (2013, 2014) and Currie et al. (2014) (~ 15 Mjup), 
this outer planet should be at least ~ 2 - 3Myr younger than 
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Fig. 12. Radial variations of the spectral index, calculated be¬ 
tween 1 and 3 mm, when a single inner planet of 20 Mjup mass 
is embedded in the disks and when two planets are in the disk, 
but the outer planet (15 Mjup) is injected 3Myr after the inner 
planet (20Mjup). The vertical lines correspond to the positions of 
the pressure maxima at the outer edge of the gap(s) carved by 
the planet(s). 

the inner planet, supporting the hypothesis that the outer planet 
may be in the act of formation. In this case, the models produce 
an outer ring of emission that is 100 times fainter than the in¬ 
ner ring of emission as observed with ALMA Cycle 0 (Fig 10). 
If the outer planet is embedded in the disk at the same time as 
the inner planet, inward drifting grains are efficiently trapped 
in the outer bump, and the contrast of the millimetre fluxes be¬ 
tween the two pressure traps becomes too strong compared with 
the observations. The ATCA observations at 3 and 7 mm also 
favour the two-planets scenario. Future high angular resolution 
and sensitivity observations of multiple transitions of optically 
thin emission from CO isotopologues, which can constrain better 
the disk gas surface density and temperature profiles, and con¬ 
tinuum images with ALMA will allow us to resolve potential 
rings (Fig. 11), to measure spectral index variations inside and 
outside pressure maxima locations (Fig. 12), and to confirm par¬ 
ticle trapping by one or two planets in this disk. 
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